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ABSTRACT 

Baryon acoustic oscillations (BAO) in the galaxy power spectrum allows us to extract the scale of 
the comoving sound horizon at recombination, a cosmological standard ruler accurately determined by 
the cosmic microwave background anisotropy data. We examine various issues important in the use of 
BAO to probe dark energy. We find that assuming a flat universe, and priors on f2 m , fl m h 2 , and Q b h 2 
as expected from the Planck mission, the constraints on dark energy parameters (wo,w') scale much 
less steeply with survey area than (area) -1 / 2 for a given redshift range. The constraints on the dark 
energy density px{z), however, do scale roughly with (area) -1 / 2 due to the strong correlation between 
H(z) and fl m (which reduces the effect of priors on O m ). Dark energy constraints from BAO are very 
sensitive to the assumed linear scale of matter clustering and the redshift accuracy of the survey. For 
a BAO survey with 0.5 < z < 2, a(R) — 0.4 (corresponding to k max {z = 0) = 0.086 h Mpc -1 ), and 
a z /(l + z) = 0.001, (a WQ ,a w >) =(0.115, 0.183) and (0.069, 0.104) for survey areas of 1000 (deg) 2 and 
10000 (deg) 2 respectively. We find that it is critical to minimize the bias in the scale estimates in order 
to derive reliable dark energy constraints. For a 1000 (10000) square degree BAO survey, a la bias in 
In H(z) leads to a 2a (3a) bias in w' . The bias in w' due to the same scale bias from IxiDa(z) is slightly 
smaller and opposite in sign. The results from this paper will be useful in assessing different proposed 
BAO surveys and guiding the design of optimal dark energy detection strategies. 
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1. INTRODUCTION 

The most intriguing mystery in cosmology today is the 
nature of the dark energy that is driving the acceler- 
ated expansion of the universe as evidenced by supernova 
(Riess et al. 1998; Perlmutter et al. 1999) and other 
data. Although current observational data are consistent 
with a cosmological constant (see for example, Wang & 
Tegmark (2004, 2005); Daly & Djorgovski (2005); Jassal, 
Bagla, & Padmanabhan (2005)), viable dark energy mod- 
els abound. See for example, Freese et al. (1987); Linde 
(1987); Peebles & Ratra (1988); Wetterich (1988); Frieman 
et al. (1995); Caldwell, Dave & Steinhardt (1998); Sahni & 
Habib (1998); Parker & Raval (1999); Dvali, Gabadadze, 
& Porrati (2000); Mersini, Bastero-Gil, & Kanti (2001); 
Freese & Lewis (2002); Carroll et al. (2004); Onemli & 
Woodard (2004); Alam & Sahni (2005); Cai, Gong, & 
Wang (2005); Cardone et al. (2005); Kolb, Matarrese, & 
Riotto (2005); Martineau & Brandenberger (2005). Pad- 
manabhan (2003) and Peebles & Ratra (2003) contain re- 
views of many models. 

The powerful complementarity of the cosmic microwave 
anisotropy (CMB) and galaxy survey data in precision cos- 
mology has long been noted (see for example, Bahcall et 
al. (1999), Eisenstein, Hu, & Tegmark (1999), and Wang, 
Spergel, & Strauss (1999)). An important development 
in this complementarity is to use the baryonic acoustic 
oscillations (BAO) in the galaxy power spectrum as a cos- 
mological standard ruler to probe dark energy (Blake & 
Glazcbrook 2003; Seo & Eisenstein 2003), made possible 
by accurate determination of the standard ruler scale by 
CMB data from WMAP (Bennett et al. 2003; Spergel et 



al. 2003) and Planck 2 . 

The "wavelength" of the baryonic acoustic oscillations 
in £;-space is determined by Ua = 2tt/s, where s is the 
comoving sound horizon at the drag epoch (when baryons 
are released from the Compton drag of the photons 3 ), 

c s dt/a = H Q ~ 1 dz—^- (1) 

10 Jz d n \ z ) 

The sound speed is 
1 

Ci- v/3(T+W 

R b = ^Pb = 31.5Q b h 2 (T CMB /2.7K)-\z/1000)-\ (2) 
and the Hubble parameter is 

h{z) = y/n m (i + zf + n rad (i + zY + n fc (i + z) 2 + n x x( z ), 

(3) 

with X{z) = px{z)/ px (0) denoting the dark energy den- 
sity function. For a cosmological constant, X(z) — 1 and 
Qx = f^A- Since Zd 3> 1, we have 

. s , / da 

y/a + a eq 
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where "d" and "eq" refer to the drag epoch and matter- 
radiation equality respectively, with Zd ~ 1089, and z eq = 
2.5 x 10 i n m h 2 (T CM B/2.7K)~ 4 . Hence the "wavelength" 
of the baryonic oscillations depends very strongly on f2 m , 
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3 For fij,h 2 < 0.03, the drag epoch follows the last scattering of the photons (Eisenstein & Hu 1998). 
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but has negligible dependence on dark energy (which only 
became important recently, at redshifts a few and less) . It 
can therefore be used as a standard ruler to probe dark 
energy. 

The systematic effects of BAO as a standard ruler are: 
bias between luminous matter and matter distributions, 
nonlinear effects, and redshift distortions (Blake & Glaze- 
brook 2003; Seo & Eisenstein 2003). Cosmological N-body 
simulations are required to quantify these effects (Angulo 
et al. 2005; Seo & Eisenstein 2005; Springel et al. 2005; 
White 2005). Fisher matrix formalism is useful in investi- 
gating the impact of various effects and assumptions that 
can be parametrized in the observed galaxy power spec- 
trum (Seo & Eisenstein 2003). 

In this paper, we examine various issues important in 
the use of baryon acoustic oscillations (BAO) to probe 
dark energy using the Fisher matrix formalism. In par- 
ticular, we extend the work of Seo & Eisenstein (2003) 
by deriving the biases in dark energy parameters (wo, w') 
due to systematic errors in the estimated scales. Through- 
out this paper, we assume spatial flatness as motivated by 
inflation and consistent with current CMB data. 

Sec. 2 describes the method used in our calculations. We 
present results in Sec. 3 and summarize in Sec. 4. 

2. THE METHOD 

2.1. Error estimation using Fisher matrix 

In estimating the expected errors in dark energy param- 
eters, we follow Seo & Eisenstein (2003) in using the Fisher 
matrix formalism (see discussion in Sec. 4) . We provide de- 
tails from our calculations for the convenience of readers 
who wish to reproduce our results. 

The comoving sizes of an object or feature at redshift 
z in line-of-sight (rii) and transverse (rj_) directions are 
related to the observed sizes Az and A9 by the Hubble 
parameter H{z) and angular diameter distance Da(z): 
rAz 

r \\ = H^-y r ± = (l + z)D A (z)A8, (5) 

The true scale of the baryonic acoustic oscillations (the co- 
moving sound horizon at recombination) is known. Hence 
if we can measure the "wavelength" of the baryonic acous- 
tic oscillations in the radial direction in successive redshift 
slices, we obtain estimates of the cosmic expansion history 
as a free function of z. While the measurement of the 
"wavelength" of the baryonic acoustic oscillations in the 
transverse direction gives us an estimate of the angular 
diameter distance D A (z) as a free function of z. 

The accuracy of the power spectrum measurement 
(Feldman, Kaiser, & Peacock 1994; Tegmark 1997) is 
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"wavelength" of the baryonic oscillations can be recovered. 
Note that in our calculation, Eq. (6) is only used to get an 
estimate of how accurately the power spectrum is deter- 
mined, and not used to derive constraints on dark energy 
and cosmological parameters. 

Assuming that the likelihood function for the band pow- 
ers of a galaxy redshift survey is Gaussian, the Fisher ma- 
trix can be approximated as (Tegmark 1997) 
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where the derivatives are evaluated at parameter values of 
the fiducial model and V e ff is the effective volume of the 
survey 
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where the comoving number density n has been taken to 
be constant in position. Here /i — k • f/fc, with f denoting 
the unit vector along the line of sight; k is the wavevec- 
tor with |k| = k. Following Blake & Glazebrook (2003) 
and Seo & Eisenstein (2003), we take k m i n = 0, and k max 
given by requiring that a(R) < 0.5 for R = n/ (2k max ) (to 
ensure that we are only considering the linear regime). 

The observed power spectrum is reconstructed using a 
particular reference cosmology, including the effects of bias 
and redshift distortions (Seo & Eisenstein 2003): 



Pobs(k re fj_, fc re /||) 
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[D A { Z )f H(z)ref 
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G(z = 0) 



.(k\z = 0) + P shot 



The values in the reference cosmology are denoted by the 
subscript "ref" , while those in the true cosmology have no 
subscript. Note that 



k re f± = k±D A (z)/D A (z) ref , k ref \ 



The 
bias 



(3 is 



= k\\H{z) ref /H{z). 

(10) 

computed from the 



where P is the average band power, V survey is the total 



rvey 

survey volume, Afc is the range of wavenumber k averaged 
over, A/j, is the range of the cosine of the angle between the 
wavevector k and the line of sight, and n is the comoving 
number density of observed galaxies. It is reasonable to 
assume nP ~ 3 (Blake & Glazebrook 2003; Seo & Eisen- 
stein 2003). This gives an estimate of how accurately the 

4 This shot noise is the unknown white shot noise that remains even after the conventional shot noise of inverse number density has been 
subtracted (Seo & Eisenstein 2003). These could arise from galaxy clustering bias even on large scales due to local bias (Scljak 2000). 



linear redshift distortion 

b for fiducial values of the observed galaxy cluster- 
ing, (3 = n^ 6 /b. G(z) is the linear growth factor. We 
normalize the power spectrum to CMB data from COBE 
(Eisenstein & Hu 1998). 

Care needs to be taken in using Eq.(9) to compute the 
derivatives of P(k) needed for the Fisher matrix in Eq.(7). 
Note that P ma tter{k\z = 0) depends on H(z) and D A (z) 
through k (see Eq.(10)). 

For a redshift slice with mean redshift z, the estimated 
parameters (assumed to be constant in the redshift slice) 
are the Hubble parameter H(z), angular diameter distance 
D A {z), linear redshift distortion (3, linear growth function 
G(z), and an unknown shot noise P s hot- 4 These are es- 
timated simultaneously with O m , U m h 2 , and Qbh 2 . The 
total number of parameters is 5N + 3 for a BAO survey 
divided into N redshift slices. 
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Note that the same cutoff scales k max (z) need to be 
used for computing all the Fisher matrix elements, in- 
cluding those of fi m , Q m h 2 , and Qi,h 2 . To implement 
this, calculate Fij in each redshift slice using the appro- 
priate k max (z), then sum over all the redshift slices for 
fl m , £l m /i 2 , and fl^h 2 . If a fixed cutoff scale is used for 
computing for the cosmological parameters, the de- 
generacy between (H(z), Da(z)) from the redshift slices 
and the cosmological parameters (especially TO ) will be 
artificially broken, leading to significant under-estimates of 
errors in (H(z), Da{z)) and the dark energy parameters. 

To obtain constraints on (w , w' , fi m , Q m h 2 ), the pa- 
rameters of interest are p = (In if 1 , h\D l Al i=l,2, N; 
fl m , fl m h 2 ). First marginalize over (/3, G(z), P s hot) in 
each redshift slice as well as VL b h 2 by taking the submatrix 
(of the parameters of interest) of the inverse of the full 
Fisher matrix, then invert it to obtain the Fisher matrix 
of the parameters of interest, F sub . The Fisher matrix 
of q = (wq, w' , f2 m , fl m h 2 ) is obtained by equating the 
log likelihood functions (ln£(p) = ln£(q)), then taking 
derivatives with respect to q on both sides. This gives 



(ii) 



2.2. 



Bias in dark energy parameters 

We now extend the work of Seo & Eisenstein (2003) by 
deriving the bias in dark energy parameters due to system- 
atic errors that bias the extracted standard ruler scale. 

The standard ruler scale is measured using H(z) and 
Da(z) in each redshift slice. First, we marginalize over 
all other parameters by taking the (In ii 4 , lnD\, i=l,2, 
N) submatrix of the inverse of the full Fisher matrix, 
then invert it to obtain the Fisher matrix relevant for scale 
determination, F scale . 

To compute the bias in (wq, w') due to biases in (In if 1 , 
hi-D^, i=l,2, N), find the (w , w') Fisher matrix by 
contracting F scale (IniP, lnD A , i=l,2, N) to 
(wo, w') (see derivation of Eq.(ll)): 

j^scate 

F ^=1^^-(F ) a0 -^- (12) 

To relate the biases in (w , w') to biases in (In if*, 
\nD A , 2=1,2, N), equate the log likelihood functions 

ln £ S ca;e( p _ p /J = l n £scale( q _ q /J + const ^ subtract 

ln£ sca ' e (p — p m ) = ln£ sca/e (q - q m ) from it, where 
Pm = Pm + s Pm, an d q' m = q m + Sq m arc the biased 
mean values. This gives 



(q - qJf Fde% (Jqjj = (p pJI (F scale ) a , VpJp 

+const. (13) 
Taking derivative with respect to qi on both sides gives 

F S DE%{^m)i = (F scale ) af} (6p m )p. (14) 



Hence 



dqi 

(HJi = (FdT);- ^ (F scale ) af3 (Spjp, (15) 

where Sp m are the biases in (lniP, \nD A , i=l,2, N), 
and summation is implied over repeated indices. 

5 This is a small effect when nP = 3 and the cosmic variance error dominates. 



A known example of BAO systematic bias is a bias in 
the dilation parameter a — k re f /ktrue for spherically aver- 
aged galaxy power spectrum. Seo & Eisenstein (2005) has 
shown that a is biased slightly above 1 in currently used 
methods for accounting the erasure of baryonic features 
due to nonlinear effects (see also White (2005)). 

3. RESULTS 

We consider a BAO survey in the redshift range of 
0.5 < z < 2, which can be carried out by obtaining 
the spectra of Ha emission line galaxies with spectro- 
graphs covering the wavelength range of 1-2 ^.m. Note that 
0.5 < z < 2 is the redshift range with the most sensitivity 
to the time variation of dark energy (Glazebrook & Blake 
2005), and easily accessible by a space mission with sim- 
ple spectroscopic instrumentation. Examples include the 
BOP MIDEX concept (Glazebrook et al. 2005), and the 
JEDI mission concept for JDEM (Wang et al. 2004; Crotts 
et al. 2005). 

Note that our results should qualitatively apply to 
other BAO surveys, in particular, ground based surveys. 
Ground based BAO surveys, for example, HETDEX (Hill 
et al. 2004) and WFMOS (Glazebrook et al. 2005), have 
been planned and will likely occur before a BAO survey 
from space. Both ground and space BAO surveys will 
be needed to establish BAO as a dark energy probe, and 
to obtain accurate and high precision dark energy con- 
straints. 

The fiducial cosmological model we have assumed is 
n s = l,h = 0.65, n m = 0.3, n A = 0.7, n b = 0.05. We as- 
sume the following priors as expected from CMB data from 
the Planck mission: <rn m = 0.01, a^ mh 2 / (fl m h 2 ) — 0.01, 
and <Jn b h2 1 (Qbh 2 ) = 0.01. These parameters are estimated 
simultaneously with H(z), Da{z), f3, G{z), and P s hot from 
seven redshift slices. Each redshift slice has the thick- 
ness Az = 0.2, except for the first redshift slice which has 
Az = 0.3 (0.5 < z < 0.8). We assume a cutoff scale k max 
given by k max = w/(2R), with cr(R) = 0.4 (except for 
Fig. 2, where a(R) is varied between 0.3 and 0.5), and a 
survey redshift accuracy of a z /(l + z) = 0.001 (except for 
Fig. 3, where a z / (1 + z) is varied between and 0.02). We 
have conservatively taken the galaxy bias to be 1 (higher 
bias leads to smaller errors in estimated dark energy pa- 
rameters 5 ), and assumed the comoving number density of 
observed galaxies such that nP = 3 following Blake & 
Glazebrook (2003). 

The dark energy equation of state is parametrized by 
w(z) = wo + w' z for z < 2. It is better to constrain dark 
energy density px{z) instead of dark energy equation of 
state (Wang & Garnavich 2001; Tegmark 2002; Wang & 
Freese 2004). This is because dark energy density is more 
closely related to observables and thus more tightly con- 
strained. We use dark energy equation of state parameters 
in this paper for easy comparison with the work by oth- 
ers. Note that the dark energy observables here are H(z), 
Da(z) from each redshift slice. The errors in H(z), Da(z) 
are propagated into the errors in Wo and w' by contracting 
the covariance matrix (see Eq.(ll)). We give constraints 
on px{z) in Fig. 5. 

Fig.l shows the la errors in wq and w' as function of 
survey area (0.5 < z < 2). In the absence of the priors, the 



4 



Yun Wang 



error bars on dark energy parameters scale as (area)" 1 / 2 
for a given redshift range. Note that with Planck priors 
on cosmological parameters, the constraints on dark en- 
ergy parameters scale much less steeply with survey area 
than (area)" 1 / 2 for a given redshift range. For a BAO 
survey with 0.5 < z < 2, (a Wo , a w i ) = (0.115,0.183) and 
(0.069,0.104) for survey areas of 1000 (deg) 2 and 10000 
(deg) 2 respectively. The smaller the survey area, the 
greater the impact of the priors on cosmological param- 
eters. 

Dark energy constraints from BAO are very sensitive to 
the linear scale of matter clustering. Fig. 2 shows the la 
errors in wq and w' as functions of a(R) for survey areas of 
1000 and 10000 (deg) 2 . The lowest panel in Fig.2 gives the 
correspondence between a(R) and the cutoff wavenumbcr 
k m ax at z = 0. Note that it is important to compare results 
for the same k max at z = 0, since the same k max (z = 0) 
may correspond to different <j(R) if different normaliza- 
tions are used for the power spectrum. 

Dark energy constraints from BAO are very sensitive to 
the redshift accuracy of the survey. Fig. 3 shows the la 
errors in w and w' as functions of a z j (1 + z) for survey 
areas of 1000 and 10000 square degrees. The redshift un- 
certainty a z leads to a radial smearing a r = ca z /H(z); 
which modifies the matter power spectrum in Eq.(9) as 
follows (see for example Peacock (1999)): 

Pobs(k\cr z ) = P obs (k\a z = 0) cxp (-fcfa 2 ) . (16) 

Fig. 4 shows an / H and ud a /Da in each of the redshift 
slices, for a z /(l + z) = 0.001, 0.01, and 0.02, for 1000 and 
10000 (deg) 2 surveys with 0.5 < z < 2. These represent 
the accuracy with which model-independent constraints 
on dark energy can be placed. These are also relevant for 
estimating likely biases in (w ,w r ) due to biases in scale 
estimates. 

Fig. 5 shows the constraints on the dark energy density 
X(z) = px{z)/ px(0) corresponding to Fig. 4 (with the 
same assumptions and the same line types). The errors 
on X(z) have been propagated from the errors on H(z) in 
each redshift slice and the error on O m (see Eq.[3]), with 
the covariance between H{z) and Q TO included. Note that 
X{z) does scale roughly as (area)" 1 / 2 , because the strong 
covariance between Q m and H(z) makes the effects of the 
cosmological priors less important. 

Table 1 gives biases in (wo,w') due to biases in (InH 1 , 
hi-D^) for 1000 and 10000 square degree BAO surveys 
with <j(R) = 0.4. Note that "lcr" denotes setting the bias 
in (lniT, lnD A ) to la errors in (lniT, \nD A ). The sign 
of the biases in (h\H l , \nD l A ) in Table 1 are chosen to cor- 
respond to scale biases in the same direction, the direction 
of wavenumber dilation a = k re f/kt rU e > 1 as found by 
Sco & Eisenstein (2005). 

If the dilation is assumed to be uniform in all directions, 
then the bias in \nH(z) and In Da(z) will be equal and op- 
posite in sign, leading to minimized biases in Wq and w' . 
However, the bias in the estimated \nH(z) should depend 
on the modeling of redshift distortions and on the redshift 
accuracy of the survey Hence it is likely that dilation in 
wavenumber will not be uniform in all directions. For this 
reason, we give the biases due to In H(z) and In Da(z) sep- 

6 Once realistic N-body simulations are available for calibrating the an 
dark energy parameters by using the data from the quasi-linear regime 



arately in Table 1. Cosmological N-body simulations are 
required to quantify the expected biases in In H(z) and 
\nD A {z). 

4. DISCUSSION AND SUMMARY 

Baryon acoustic oscillations (BAO) provide an impor- 
tant new method for probing dark energy To extract reli- 
able dark energy constraints from BAO surveys, we need 
to understand how the expected errors depend quantita- 
tively on the assumptions made. 

For simplicity and transparency, we use the Fisher in- 
formation matrix formalism in this paper. Our results are 
consistent with those of Seo & Eisenstein (2003), and agree 
to about 10% with those found using Monte Carlo methods 
by Blake & Glazebrook (2003) and Glazcbrook & Blake 
(2005), for the same k max (z — 0) and priors on cosmolog- 
ical parameters. Although Fisher matrix yields the small- 
est possible error bars, we have made the most conservative 
assumptions. In deriving the dark energy constraints, we 
marginalize over many additional physical parameters (Seo 
& Eisenstein 2003): the linear redshift distortion (5, linear 
growth function G(z), and an unknown shot noise P s hot 
in each redshift slice; the matter density ratio fi m , the 
matter density Q, m h 2 , and the baryon density flbh 2 . This 
approach seems to work rather well, as Seo & Eisenstein 
(2005) found that their results from cosmological N-body 
simulations are close to what they found using the Fisher 
matrix formalism (Sco & Eisenstein 2003). 

We have examined various issues important in the use 
of BAO to probe dark energy. In particular, we extend the 
work of Seo & Eisenstein (2003) by deriving the biases in 
dark energy parameters (w , w') due to systematic errors 
in the estimated scales. 

We find that assuming priors on Ct m , fl m h 2 , and flbh 2 as 
expected from the Planck mission, the constraints on dark 
energy parameters (w ,w r ) scale much less steeply with 
survey area (for a given redshift range) than (area)" 1 / 2 
(which holds in the absence of the priors). For a BAO 
survey with 0.5 < z < 2, a(R) = 0.4 (corresponding to 
k m ax(z = 0) = 0.086 hMpc- 1 ), and a z /(l + z) = 0.001, 
{(T W0 ,a w >) = (0.115,0.183) and (0.069,0.104) for survey 
areas of 1000 (deg) 2 and 10000 (deg) 2 respectively (see 
Fig.l). This provides a useful guide in designing optimal 
survey strategies that provide robust and accurate dark 
energy constraints. 

The constraints on the dark energy density px(z), how- 
ever, do scale roughly with (area) -1 / 2 due to the strong 
correlation between H(z) and O m (which reduces the effect 
of priors on fl m ). This is interesting since px(z) provides 
direct and model-independent constraints on dark energy 
(Wang & Garnavich 2001; Tegmark 2002; Wang & Freese 
2004). 

The BAO constraints on the dark energy equation of 
state parameters (too, to') are very sensitive to the assumed 
linear scale of matter clustering k max (z) 6 and the redshift 
accuracy of the survey cr z /(l + z) (see Figs. 2-3). This is 
important to note since different proposed BAO surveys 
are often compared by their expected errors in (w ,w'). 
Such comparisons are not appropriate unless the same 
assumptions are made about k max (z), and the claimed 

sis of real data, we will be able to extract additional information on 
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a z /(l + z) can be demonstrated to be feasible. The lat- 
ter is a key issue for surveys using photometric redshifts 
(Blake & Glazebrook 2003; Seo & Eisenstein 2003; Zhan 
& Knox 2005). 

It should be noted that the radial BAO are not observ- 
able for photometric rcdshift accuracies. The H(z) mea- 
surement from using photometric redshifts is derived from 
the degree of damping of broad-band power in the radial 
direction on very large scales (k < 0.05hMpc _1 ). Hence 
this measurement is much less robust (more susceptible to 
systematic error) than the spectroscopic case, where the 
radial BAO are detectable and give H(z) measurements 
directly. 

We find that it is critical to minimize the bias in the 
scale estimates (lni?(z) and In Da(z)) in order to derive 
reliable dark energy constraints. For a 1000 (10000) sq 
deg BAO survey, a la bias in ln_ff(z) leads to a 2a (3a) 
bias in w' . The bias in w' due to the same scale bias 
from ItiDa(z) is slightly smaller and opposite in sign. In 
addition to the scale bias due to the inaccurate model- 
ing of nonlinear effects (which impacts the results in the 
linear regime, sec Seo & Eisenstein (2005)), systematic 
biases in the H(z) measurement could arise when photo- 



metric redshifts are used. It is important to quantify the 
expected bias in h\H(z) and \t\Da{z) using cosmologi- 
cal N-body simulations. Note that H(z) and Da{z) are 
measured as independent parameters. Since Da{z) is re- 
lated to H (z) through an integral for a given cosmological 
model in the absence of systematic biases, comparing the 
directly measured Da(z) to the Da{z) derived from the 
H{z) measurement provides us with a cross-check to help 
model unknown systematic effects. 

The planned BAO surveys from both ground and space 
telescopes will play an important role in unraveling the 
nature of dark energy. The results from in this paper will 
be useful in assessing different proposed BAO surveys and 
guiding the design of optimal dark energy detection strate- 
gies. 
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Table 1 

Biases in (wq,w') due to biases in (lnJT, lnD\). 



survey area 


bias in In Da{z) 


bias in In H(z) 


o z /{l + z) 


bias in wq (a Wo ) 


bias in w' (a w /) 


1000 (deg)^ 


-la 


0. 


0.001 
0.01 
0.02 


0.082 (0.115) 
0.255 (0.266) 
0.273 (0.396) 


-0.375 (0.183) 
-0.930 (0.727) 
-1.021 (1.169) 




0. 


la 


0.001 
0.01 
0.02 


-0.181 (0.115) 
-0.303 (0.266) 
-0.421 (0.396) 


0.317 (0.183) 
0.491 (0.727) 
0.742 (1.169) 


10000 (deg) z 


-la 


0. 


0.001 

0.01 

0.02 


0.033 (0.069) 
0.085 (0.133) 
0.089 (0.167) 


-0.291 (0.104) 
-0.311 (0.251) 
-0.317 (0.386) 




0. 


la 


0.001 
0.01 
0.02 


-0.078 (0.069) 
-0.137 (0.133) 
-0.157 (0.167) 


0.263 (0.104) 
0.172 (0.251) 
0.198 (0.386) 



Dark Energy and Baryon Acoustic Oscillations 
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Fig. 1. — The la errors in wq and w' as function of survey area (0.5 < z < 2). Note that the constraints on dark energy parameters 
much less steeply with survey area (for a given redshift range) than (area) -1 ' 2 (which holds in the absence of the priors). 
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Fig. 2. — The 1<t errors in ivq and w' as function of &(R) for survey areas of 1000 and 10000 square degrees. The lowest panel gives the 
correspondence between cr(ij) and the cutoff wavenumber k ma x at Z = 0. 
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Fig. 3. — The lcr errors in wq and to' as function of cr z /(l + z) for survey areas of 1000 and 10000 square degrees. 
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Fig. 4. — The estimated cru/H and od a /Da in each of the redshift slices, for cr z /(l + z) = 0.001, 0.01, and 0.02 for survey areas of 1000 
and 10000 square degrees. 
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Fig. 5. — The estimated dark energy density X(z) = px(z)/px (0) corresponding to Fig. 4 (with the same assumptions and the same line 
types). 



